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We calculate the thermal conductivity of electrons produced by electron-electron Coulomb scat- 
tering in a strongly degenerate electron gas taking into account the Landau damping of transverse 
plasmons. The Landau damping strongly reduces this conductivity in the domain of ultrarelativistic 
electrons at temperatures below the electron plasma temperature. In the inner crust of a neutron 
star at temperatures T < 10 7 K this thermal conductivity completely dominates over the electron 
conductivity due to electron-ion (electron-phonon) scattering and becomes competitive with the the 
electron conductivity due to scattering of electrons by impurity ions. 
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I. INTRODUCTION 

The electron thermal conductivity is an important ki- 
netic property of plasmas in experimental devices, in 
metals and semiconductors and in astrophysical objects. 
It has been studied for a long time and described in text- 
books (see, e.g., Ziman 0], Lifshitz and Pitaevskii 0]). 
We will show that some aspects of this problem have to 
be revised. 

Specifically, we will consider the thermal conductiv- 
ity of degenerate electrons. It is important in metals, 
semiconductors as well as in degenerate cores of evolved 
stars (giants, supergiants and, most importantly, white 
dwarfs), and in envelopes of neutron stars. It is needed 
to study cooling of white dwarfs (e.g., Prada Moroni and 
Straniero |3|) and nuclear explosions of massive white 
dwarfs as type la supernovae (e.g., Baraffe, Heger and 
Woosley Q). In neutron star envelopes it is required for 
many reasons — for simulating cooling of isolated neu- 
tron stars (see, e.g., Lattimer et al. Gnedin, Yakovlev 
and Potekhin 6], Page, Geppert and Weber @); for 
studying heat propagation and thermal relaxation in a 
neutron star envelope in response to the various processes 
of crustal energy release. In particular, these processes 
include pulsar glitches (e.g., Larson and Link and ref- 
erences therein), deep crustal heating of accreting neu- 
tron stars in soft X-ray transients (e.g., Ushomirsky and 
Rutledge 0), superbursts as powerful nuclear explosions 
at the base of the outer crust of an accreting neutron star 
je.g ., Strohmayer and Bildsten |lfj|: Page and Cumming 

It is well known, that the electrons give the main con- 
tribution into the thermal conductivity of strongly degen- 
erate matter. Their thermal conductivity can be written 
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Here, T is the temperature, fee is the Boltzmann con- 
stant, n e is the number density of electrons, m* = /i/c 2 , 
ji is the electron chemical potential (including the rest- 
mass term) , and v e is the total effective electron collision 



frequency. The latter frequency is the sum of the par- 
tial collision frequencies v e i and u ee . In v e i we include 
all collisions of electrons mediated by their interactions 
with ions (direct Coulomb scattering of electrons by ions 
in an ion gas or liquid; electron-phonon scattering in an 
ion crystal and electron scattering by impurity ions at 
low temperatures). Evidently, Eq. can be rewritten 



as 
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where the partial conductivities K e i and deter- 
mined, respectively, by the ei and ee collisions. 

It is widely believed that the dominant contribution 
into K e comes from the ei collisions. The associated 
partial conductivity K e i has been studied in a number 
of papers (e.g., Potekhin et al. 12], Gnedin et al. [f| 
and references therein). The partial thermal conductiv- 
ity K ee owing to the ee collisions was calculated by Lampe 
[l3| , Flowers and Itoh Q , Urpin and Yakovlev , and 
Timmes [T^ . These calculations were summarized by 
Potekhin et al. 01 ■ The main result was that the ee col- 
lisions are mainly negligible, except for a hot low-density 
plasma of light elements (from hydrogen to carbon) with 
the temperature a few times lower than the electron de- 
generacy temperature (in that case the ee collisions did 
affect K e but never dominated as long as the electrons 
were degenerate; see Lampe 0|; Urpin and Yakovlev 
0). 

However, all these calculations have neglected an im- 
portant effect of the Landau damping of the ee interac- 
tion owing to the exchange of transverse plasmons. In the 
context of transport properties of dense matter this effect 
was studied by Heiselberg and Pethick 0] for degener- 
ate quark plasmas. Similar effects have recently been 
analyzed by Jaikumar, Gale and Page |l8j for neutrino 
bremsstrahlung radiation in the ee collisions. Here we 
reconsider n ee including the effects of the Landau damp- 
ing and show that K ee is actually much more important 
than thought before. 
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II. FORMALISM 

We analyze the thermal conductivity n ee of strongly 
degenerate electrons under the conditions that the elec- 
trons constitute an almost ideal and uniform Fermi gas 
and collide between themselves and with plasma ions. 
Although the results will be more important for ultrarel- 
ativistic electrons, our consideration will be general and 
valid for both, relativistic and non-relativistic, cases. We 
will use the same standard variational approach with the 
simplest trial non-equilibrium electron distribution func- 
tion as Flowers and Itoh [14j . The calculations are similar 
to those performed by Heiselberg and Pethick ^3 f° r the 
thermal conductivity of quarks; thus we omit the details. 

Following Heiselberg and Pethick [l7| we have 



(i) 
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where p e = h (37r 2 rt e ) 1 / 3 and v e = p e /m* are, respec- 
tively, the electron Fermi momentum and Fermi veloc- 
ity; the integration is over all allowable ee collisions 
P1P2 —> PiP 2 ] P is an electron momentum, v its ve- 
locity, e its energy; primes refer to electrons after a col- 
lision event; / is the electron Fermi-Dirac distribution, 
and VF(12|1'2') is the differential transition probability 

VF(12|1'2') = 4(27Th) 6 S{£' 1 +£ , 2 -£ 1 ~£ 2 ) 

x5(p' l +p l 2 -p 1 ~p 2 )\M fl \ 2 , (4) 

being the squared matrix element summed over 
electron spin states. The symmetry factors required to 
avoid double counting of the same initial and final elec- 
tron states are also included into \Mfi\ 2 . 

The multi-dimensional integral (0) is simplified fur- 
ther using the standard angular-energy decomposition 
(separating integrations over particle energies and ori- 
entations of their momenta; e.g., Shapiro and Teukol- 
sky |l9j). Because the electrons are strongly degener- 
ate, they can participate in thermal conduction only if 
their energies are close to the Fermi level. Accordingly, 
their momenta can be placed at the Fermi surface in an- 
gular integrals whenever possible. Characteristic energy 
transfers ftoj = e[ — £\ in collisions of strongly degener- 
ate particles are small, huj < fceT. Momentum transfers 
hq = p[ — pi are also small, hq <C p e , owing to a long- 
range nature of the Coulomb interaction. We will use this 
small-momentum-transfer approximation throughout the 
paper. Typical values of Hq are determined by plasma 
screening of the Coulomb interaction. 

The plasma screening was thoroughly analyzed by 
Heiselberg and Pethick ^3- These authors studied 
quark-quark interaction through one-gluon exchange in 
the weak-coupling limit which is very similar to the 
Coulomb interaction in an ordinary plasma. The matrix 



element for an ee scattering event is Mfi = 
r(i) 



■M 



r (2) 

fi 

2'; 2 



(2) 
fi ' 



where M^' and M^' correspond to the channels 1 
1': 2^2' and 1 



1', respectively. For instance, 
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where 



We'e ' 



J e 'e) = c(u e '"f v u e ) is the transition 
4-current (^=0, 1, 2, 3), Jte'e is the component of J e / e 
transverse to q, j u is a Dirac matrix, u e is a normal- 
ized electron bispinor (u e u e = 2m e c 2 ), and u e is a Dirac 
conjugate (see, e.g., Berestetskii, Lifshitz and Pitaevskii 
ppf). The longitudinal component of J e < e (parallel to 
q) is related to the time-like (charge density) component 

Jg/g via current continuity equation; it is excluded from 
Eq. ©, see Ref. [13 . 

The polarization functions B; and B t in Eq. (|SJ) de- 
pend on u> and q and describe the plasma screening of 
the ee interaction through the exchange of longitudinal 
and transverse plasmons, respectively. In the classical 
limit (hq <C p e and hw <C v e p e ), taking into account 
the polarization produced by degenerate electrons in the 
first-order random phase approximation, one has (e.g., 
Alexandrov, Bogdankevich and Rukhadze |2l|]) 



n ; = ql Xl( x )> n t = [qaVe/c) 2 Xt(x), 



(6) 



where x — Lu/(qv e ), 



with 
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qo being the ordinary Thomas-Fermi electron screening 
wavenumber. Particularly, in the limit of u> — > and 
Lo/q <C v e we have 



Xi 



1, 



Xt = inuj/(4qv e ). 



(9) 



According to Eqs. (JHJ and ©, the plasma screening of 
the ee current interaction via the exchange of transverse 
plasmons (the second term in Eq. JSJ) is different from 
the screening of the charge interaction via the exchange of 
longitudinal plasmons (the first term). The difference re- 
sults from the difference of the polarization functions Tit 
and B; and has been neglected in all previous calculations 
of the electron thermal conductivity (where one has com- 
monly set B f = Bf = ql, and q 2 - (uj/c) 2 + U t = q 2 +q 2 )- 
Naturally, this difference is expected to be small for non- 
relativistic electrons (v e <C c), where the transverse cur- 
rent interaction term is small in the matrix element. We 
will see that the difference becomes important for rela- 
tivistic electrons (v e ~ c). 
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Let us recall that the functions 11/ and Il t have real 
parts which describe plasmon refraction, and imaginary 
parts which describe plasmon absorption. In the ee 
scattering we deal with low-energy virtual plasmons, 
hui < fceT. As seen from Eq. JDJ, longitudinal plas- 
mons undergo refraction which results in the Debye-type 
(Thomas-Fermi) screening of the Coulomb interaction, 
with the screening momentum (inverse screening length) 
qq. As for transverse plasmons, they mainly undergo col- 
lisionless absorption (that is the Landau damping) by 
degenerate electrons. Their effect is drastically different 
from the effect of longitudinal plasmons. 

The calculations similar to those in Ref. |l7j lead to 
the following expressions for the ee collision frequency 
and thermal conductivity, 



36n e a 2 h 2 cl(u, ( 
-Kratk^T 
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Here, a = e 2 /Tic is the fine structure constant, and 
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1 — u 2 x 2 + u 2 (x6/w) 2 Xt(x) 
is a dimcnsionless function of two variables, 



(11) 



v e /c, 



hv e q /(k B T) = \/3T pe /T, (12) 



T pe = hiu pe /kB being the electron plasma temperature 
determined by the electron plasma frequency uo pe — 
\J A'Ke 2 n e /m*. Both components of the matrix element, 

and M**', give equal contributions into Eq. (fTT|l . 
and the interference term is negligibly small because of 
the small-momentum-transfer approximation. Further- 
more, w — hui/ksT; 4> is the angle between p\ t and p2t, 
the components of pi and P2 transverse to q; and the 
integration over <f> is trivial. Equations (|10|l and (|llfl are 
natural generalizations of Eqs. (58) and (59) of Heisel- 
berg and Pethick 01 to the case of a degenerate gas 
of particles of arbitrary degree of relativity (in our case 
v/c may be arbitrary while in Ref. [T^ v — c). The dif- 
ference of numerical factors in the expressions for k in 
our Eq. and in Eq. (58) of Ref. Qj (108 versus 24) 
stems from the difference of physical systems under con- 
sideration (an electron gas versus a gas of light quarks 
interacting through gluon exchange). 



III. FOUR REGIMES OF 
ELECTRON-ELECTRON COLLISIONS 

Thus the collision frequency and the thermal conduc- 
tivity IjlOfl are solely determined by the electron number 



density and the temperature. Their calculation reduces 
to the calculation of the function I(u,9) from Eq. 111( 1 . 
Clearly, the function can be written as 



I = h+I t + In 



(13) 



where // is the contribution from the ee interaction via 
the exchange of longitudinal plasmons (the first term in 
the squared modulus); It comes from the interaction via 
the exchange of transverse plasmons (the second term); 
and In is the mixed term. 

The analysis reveals four regimes (I-IV) of ee collisions 
in a strongly degenerate electron gas. These regimes are 
summarized in Table [I] The regimes I and II are realized 
for non-relativistic electrons, while in the regimes III and 
IV electrons are ultrarelativistic. The regimes I and III 
take place for sufficiently high temperatures T > T pe , at 
which the Pauli principle does not restrict energy trans- 
fers between colliding electrons (Huj < k^T; see, e.g., 
Lampe 0], especially his Fig. 1). The regimes II and 
IV refer to a colder electron gas, where energy transfers 
are essentially limited by the Pauli blocking (e.g., Lampe 
PI and Flowers and Itoh 0). 

The analysis of Eq. I|ll|) gives the following asymptotic 
values of / in the different regimes. 

In the regime I (where u < 1 and < 1) 



h 
It 
Ili 



I (A Inl +0.1657), 
U " (315 ln ^ + °-° 5067 )' 



^1,1 +0.1236) 



(14) 



The logarithmic terms in brackets represent Coulomb 
logarithms, while the second terms are the corrections 
calculated using the standard technique The lead- 
ing contribution comes from Ii. It was calculated by 
Lampe [l3| , with a slightly less accurate correction term 
(1.30 x 2/15 w 0.173 from his Eqs. (5.22) and (5.23), in- 
stead of our 0.1657). Retaining this term in the regime 
I, one has I = Ii and 
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72a 2 ft 2 c 2 [ln(l/0) + 1.242]' 



(15) 



Table I: Four regimes of thermal conduction of degenerate 
electrons owing to ee collisions. 





Electron 




Main 


T-dependence 


Regime 


velocity 


Temperature 


contribution 


Of K ee 


I 


We < c 


T > T pe 


h 


T 2 /\n(T/T P e) 


II 


V e < C 


T< T pe 


h 


1/T 


III 


V e W C 


T > Tpe 


h + It + lit 


T 2 /\n(T/T P e) 


IV 


V e « C 


T <C T p e 


It 


const 



4 



Note that for I t the regime I extends to lower tempera- 
tures T ~ uT pe , than for and In, but this circumstance 
does not affect noticeably the thermal conductivity n ee 
because It is relatively insignificant in the given regime. 
In the regime III (where u « 1 and 9 < 1), 



I, 

hi 

I 



2/ * = 77 mi + 0.1657, 
15 t) 

^ In i + 0.1399, 
15 t; 

- In - + 0.3884. 
3 



(16) 



Again, the logarithmic terms are Coulomb logarithms 
and the second terms are corrections. In this case, all 
terms (// , It , and Iu ) give comparable contributions into 
/. The asymptote of / in this regime was obtained 
by Heiselberg and Pethick 01 (their Eq. (60), with a 
slightly less accurate correction factor 0.30 instead of our 
0.3884) in their studies of quark plasma. For an elec- 
tron plasma, the conductivity K ee in the regime III was 
calculated by Urpin and Yakovlev [l^. Their result is 
equivalent to Eq. (|16fl for / but with less accurate cor- 
rection (ln(2)/3 w 0.231 instead of 0.3884) because they 
erroneously used the approximation of static longitudi- 
nal electron screening of the Coulomb interaction in all 
terms (with II; — H t — q 2 , in Eq. ©). Employing / from 
Eq. O in the regime III we obtain 
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In the regimes II and IV (where 9 > 1 and u < 1 is 
arbitrary) we have 



h = 



" '* = 2C(3)Jr, Iu = Jr^ 1/3 , (18) 



15u6> 3 



18.52, 



where £(z) is the Riemann zeta function (with £(3) = 
1.202) and T(z) is the gamma function. The expres- 
sion for It ( with u = 1) was obtained by Heiselberg and 
Pethick |l7J for an ultrarelativistic quark plasma. 

The asymptotic expressions (|18t are derived from 
Eq. (|ll(l with the screening functions xi an d Xt given 
by Eq. (0. In this approximation the screening of the 
ee interaction via the exchange of longitudinal plasmons 
(the first part of the matrix element in Eq. JSJ) is de- 
scribed by q 2 + Hi w q 2 + q 2 7 which is equivalent to the 
static Debye-type screening with the screening wavenum- 
ber q scl -.i = qo. The screening via the exchange of trans- 
verse plasmons (the second part of the matrix element) 
is more complicated. It is described by the denominator 
term q 2 — uj 2 /c 2 +H t ~ q 2 + iq 2 nLov e / (4gc 2 ) that repre- 
sents the dynamical screening via the Landau damping 
of transverse plasmons, with hu; ~ k^T <C v e hq (i.e., 
with low phase velocities u/q <C v e ). In this case the ef- 
fective screening wavenumber q SCIt t is evidently given by 
3scr,t ~ <loUv e /c 2 ~ q 2 k B Tv e / (he 2 ). 



Using Eqs. i|10|) and (|13fl we can decompose the col- 
lision frequency v ee into the same parts as /, v ee = 
Vee + Vee ■ One can easily see, that in the regimes 
II and IV the partial frequency vie can be estimated as 
n e a 2 (k^T) 2 / (hql cr ; /i). The collision frequency 



vi} can be estimated by replacing 

9scr,z — * 9scr,t- The factor u 4 takes into account the re- 
duction of the ee interaction via transverse plasmons in 
the nonrelativistic electron gas. The quantity v£$ can 



be estimated as v, 



(») 



u 2 vi} with ql 



3 2 
iscr,l ~^ 9scr,i <Zscr,i- 

In the regime II, where 8 > 1 and u -C 1, the exchange 
of longitudinal plasmons dominates, and we have / sa //, 



bhqlv 4 
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This leading part of n ee in the regime II was calculated 
by Lampe [pj ■ Note that the asymptotic expression ijT%|) 
for If at u <C 1 is actually valid not at T < T pe , as the 
expressions for // and In, but at lower T < uT pe . How- 
ever, this circumstance is relatively unimportant because 
it is which dominates in the regime II. 

In the regime IV, where > 1 and u « 1, the exchange 
of transverse plasmons dominates, with I w I t , and 



7TfcBCgo 

216C(3)a 2 ' 



(20) 



It is remarkable that in this regime K ee becomes temper- 
ature independent. This regime has been discussed in 
detail by Heiselberg and Pethick 0] for quark plasma. 
For the electron plasma, it was considered by Flow- 
ers and Itoh [l4| and also by Urpin and Yakovlev [l^ 
but both groups erroneously used the approximation of 
static longitudinal screening in all channels (xi = Xt — 
1) which strongly underestimates the efficiency of the 
plasma screening of the ee interaction in the ultrarela- 
tivistic electron gas. If that approximation were true, 
one would obtain v, 



(0 



,(«) 



whereas in fact VeJ *C vie <C VeV in the regime IV, and 
K e e is significantly lower than predicted by the previous 
calculations [l4lll5| . 

Let us note that the temperature behavior of n ee (Ta- 
ble [U corresponds to an ordinaryFermi-liquid (where 
k cx 1/T; e.g., Baym and Pethick |22j) only in the regime 
II. In the regimes I and III the plasma is too warm (al- 
though degenerate) to reach the Fermi-liquid limit, where 
energy transfers between colliding particles are strongly 
restricted by the Pauli principle. In the regime IV the 
plasma is cold but the matrix element of the ee interac- 
tion essentially depends on energy transfers ftiu owing to 
the Landau damping of transverse plasmons which vio- 
lates the Fermi- liquid behavior [l7| . 

To facilitate applications, we have calculated the func- 
tions Ii , It and In from Eq. (|llfl for a dense grid of u and 
9 and fitted the numerical data by analytic functions. We 



1v£} and v & 
,(*) 



2.5z4i ) , 
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obtain 



- I 0.1587 

u 

x In ( 1 + 



0.02538 
1 + 0.0435 ( 
128.56 



37.1 0+ 10.83 9 2 + 9 3 



(21) 



Therefore, // is inversely proportional to u, with the pro- 
portionality coefficient dependent solely on 9. Further- 
more, 



It = u 3 



2.404 C 2 - 2.404/C 

C M 



x In 1 



1 + 0.1 6u 
C 



A9u + 9 2 



(22) 



where A = 20 + 450u 3 , C = Aexp(C x /C 2 ), C x = 
0.05067 + 0.03216 m 2 , and C 2 = 0.0254 + 0.04127 u 4 . Fi- 
nally, 



lit = u 



18.52 u 2 C 2 - 18.52 u 2 /C 



C 

x In I 1 + 



1 + 0. 1558 s 
C 



A6 + 10.83 9 2 u 2 + (9uf/ 3 



, (23) 



where A = 12.2 + 25.2 it 3 , B = 1 - 0.75 u, C 
A exp(Ci/C 2 ), Cx = 0.123636 + 0.016234 u 2 , and C 2 = 
0.0762 + 0.05714 it 4 . These fit expressions reproduce also 
all asymptotic limits mentioned above. The maximum 
fit errors of I t , and J/ t are, respectively, 0.2% (at 
u = 0.61 and 9 = 3.8), 6.3% (at u = 0.878 and 9 ~ 0.175, 
where the electrons become nondegenerate) and 8.4% (at 
u = 0.19 and 9 = 19). The maximum fit error of the to- 
tal function I in Eq. ((El is 2.5% (at u = 0.924 and 
9 = 0.232). 



IV. DISCUSSION 

Let us discuss the efficiency of ee collisions for thermal 
conduction in a degenerate electron gas. Figure ^ shows 
the temperature-density diagram for a dense matter. For 
certainty, we have adopted the model of cold-catalyzed 
(ground-state) matter. It is composed of electrons and 
ions (atomic nuclei) at densities smaller than the neutron 
drip density p n d ~ 4x 10 11 g cm -3 (the right vertical dot- 
ted line), with an addition of free neutrons at higher den- 
sities p. The composition of the cold-catalyzed matter is 
taken from Haensel and Pichon j2^| at p < p n d and from 
Negele and Vautherin [24j at higher p. Weak first-oder 
phase transitions which accompany changes of ground- 
state nuclides with the growth of p are smoothed out as 
described by Kaminker et al. [2{|. At p < 10 8 g cm~ 3 
the ground-state matter is composed of iron ( 56 Fe) . The 
densities p < 10 9 g cm -3 are appropriate for degenerate 
stellar cores of white dwarfs and evolved stars; the den- 
sities p < 10 14 are appropriate for envelopes (crusts) of 
neutron stars. 
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Figure 1: (color online) Temperature-density diagram for 
cold-catalyzed matter. We show the electron degeneracy tem- 
perature Tp, the electron plasma temperature T pe , the melt- 
ing temperature T m of the ion crystal, and the ion plasma 
temperature T pi . The left dotted line separates the regions 
of a low-density non-relativistic electron gas and a denser gas 
of ultrarelativistic electrons. The right dotted line indicates 
the neutron drip point. Shaded regions I-IV show the T — p 
domains of the different ee collision regimes (Tabic QJ. 



In Figure we plot the electron degeneracy temper- 
ature Tp = (p — m e c 2 )/k-Q, the electron plasma tem- 
perature T pe , the melting temperature of the ion crystal 
T m ~ Z 2 e 2 1 '(175afcB), and the ion plasma temperature 
Tpi = Hu>pi/kB (where io V i — yj 4TrZ 2 e 2 rii /rrii is the ion 
plasma frequency, Zie is the ion charge, m, is its mass, 
n.j is the ion number density, and a = (47rnj/3) -1 / 3 is 
the ion-sphere radius) . The left vertical dotted line sepa- 
rates the regions of non-relativistic degenerate electrons 



( V e 



C, p > 



>6 

10 6 



(v e <C c, p <C 10° g cm 3 ) and ultrarelativistic electrons 



g cm 3 ). The shaded regions I-IV 



show the temperature-density domains with the different 
regimes of ee collisions in a degenerate electron gas (see 
Section ITTT1 Tabled). 

In order to explore the importance of the ee collisions 
for thermal conduction let us compare K ee with the ther- 
mal conductivity n e i owing to ei collisions (Section QJ- 
The conductivity K e i is calculated using the formalism of 
Potekhin et al. jTj and Gnedin et al. |||. 

In Figure|21we plot the temperature dependence of the 
electron thermal conductivity of helium ( 4 He) or carbon 
( 12 C) plasmas at p = 10 5 g cm -3 (left) or 10 6 g cm~ 3 
(right). At p = 10 5 g cm" 3 the electrons are nonrelativis- 
tic while at p = 10 6 g cm -3 they are mildly relativistic. 
According to Eq. , the total electron thermal conduc- 
tivity n e is determined by the partial contributions K ee 
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Figure 2: (color online) Temperature dependence of the elec- 
tron thermal conductivity at p = 10 5 g cm -3 (left) and 10 6 
g cm -3 (right). The dashed line ee shows « ee ; the dotted line 
ee — Hs the same but retaining the contribution of longitudi- 
nal plasmons alone. The dot-and-dashed lines show n e i and 
the solid lines show the total (tot, ee + ei) conductivity re e for 
helium and carbon plasmas. 



and K e i, a minimum partial contribution being most im- 
portant. The total conductivities are shown by the solid 
lines; K e t by the dot-and-dashed lines; n ee is plotted by 
the dashed lines (the same for the helium and carbon 
plasmas); the dotted lines give K ee neglecting the contri- 
bution of transverse plasmons (i.e., by setting I = 

For the densities and temperatures displayed in Figure 
121 the ei collisions are mainly more efficient than the ee 
ones, although the contribution of the ee collisions into 
the total thermal conductivity is noticeable. One can 
see that the ei collisions are more important for heavier 
elements (because the Coulomb ei scattering cross sec- 
tion is much higher than the ee scattering cross section 
for high-i? elements; see Lampe 01). The ee collisions 
would be negligible for the iron plasma if the data for this 
plasma were displayed in Figure |21 The ee collisions are 
more efficient in the helium plasma than in the carbon 
plasma because the helium ions have lower charges. For 
a given chemical composition, n ee gives highest contri- 
bution into K e at temperatures T a few times lower than 
T pe (we have log 10 T pe [K] w 7.86 and 8.32 at p = 10 5 and 
10 6 g cm -3 , respectively). These temperatures separate 
the high-temperature and low-temperature ee collision 
regimes (e.g., the regimes I and II in the non-relativistic 
electron gas, see Table QJ. For p = 10 5 g cm~ 3 , the elec- 
tron gas is only slightly relativistic. Accordingly, the con- 
tribution of the Landau damping (transverse plasmons) 
into K ee is relatively small, and the results of Lampe [l3| 
are sufficiently accurate. For higher p = 10 6 g cm -3 , 
the contribution of the Landau damping becomes more 
important which invalidates the previous results of Refs. 
[l4L fl5l | . In this case the Landau damping strongly re- 
duces K ee which makes the ee collisions much more im- 
portant for thermal conduction than it was thought be- 
fore. For still higher p the contribution of n ee into K e 
in a plasma of light elements at T ~ T pe would be even 
more significant. However, the temperature T pe would 
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Figure 3: (color online) Density dependence of the electron 
thermal conductivity at different T (marked by the values of 
log 10 T near the curves) in a neutron star crust composed of 
ground-state matter; the vertical dotted line is the neutron 
drip point. Dot-and-dashed lines (ei) show the conductivity 
K e i owing to ei collisions at T = 10 7 , 10 s , and 10 9 K in 
a crystal of atomic nuclei for pure ground-state matter. In 
addition, the lower dot-and-dashed line at T = 10 7 K shows 
K e i which includes the contribution of electron scattering by 
impurity nuclei at T — 10 7 K (for ni mp /ni = 0.05 and |Zi mp — 
Zi\ =2). The dashed line is K ee ; it is almost independent of 

T. The dotted line is the partial conductivity nil mediated 
by longitudinal plasmons alone at T — 10 8 K; it scales as 
T~ . The solid lines show the total electron conductivity K e 
at T — 10 7 K for the matter without and with impurities. If 
T > 10 8 K the effect of impurities is weak and K e ~ K e %- 



become so high that light elements would start burning 
in thermonuclear reactions. Therefore, ee collisions can 
noticeably decrease the electron thermal conductivity at 
densities and temperatures important for nuclear burn- 
ing of light elements (for instance, in the vicinity of the 
carbon ignition curve). This should be taken into ac- 
count in simulations of nuclear evolution of stars (e.g., 
carbon ignition in white dwarfs Q). 

Figure |31 shows the density dependence of the elec- 
tron thermal conductivity in the density range from 
~ 2.5 x 10 9 g cm -3 to 10 14 g cm -3 for the three val- 
ues of the temperature T = 10 7 , 10 8 , and 10 9 K. We 
employ the same model of ground-state matter as in Fig- 
ure n The displayed density range is appropriate to a 
crust of a neutron star. The vertical dotted line is the 
neutron drip point which separates the crust into the 
outer and inner parts. The adopted values of T and p 
refer to the regime IV of the ee collisions, where the elec- 
tron gas is ultrarelativistic, T <C T pe , and the Landau 
damping dominates. The conductivity K ee (the dashed 
line) is well approximated by the temperature indepen- 
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dent conductivity (|20() governed by the Landau damping; 
K ee slightly deviates from Eq. I|20|) only at lowest p for 
T = 10 9 in Figure where T is still insufficiently lower 
than T pe (see Figure Retaining the longitudinal con- 
tribution into K ee (I = Ii) and setting T = 10 8 K, we 
would get the dotted curve in Fig. ft goes much higher 
than the dashed curve indicating that this longitudinal 
contribution is really insignificant. 

For comparison, the dot-and-dashed lines in Figure |3 
give the thermal conductivity K e i (mostly from Figure 4 
of Ref. 6]) for the same values of T, and the solid lines 
give the total conductivity K e . For T = 10 s and 10 9 , 
the effects of possible impurities in dense matter (atomic 
nuclei with charge numbers Z- lmp different from charge 
numbers Zi of ground-state nuclei) are expected to be 
small (e.g., Gnedin et al. @). At T = 10 7 K, the effects 
of impurities can be substantial. For this T in Figure 
OH we present n e i and K e for pure ground-state matter 
and also for the matter which contains impurities with 
l^imp — Zi\ = 2 and with the fractional number of im- 
purity nuclei of ni mp /rii = 0.05. The impurities increase 
the ei collision rate and decrease the conductivities K e j 
and n e . 

As seen from Figure |3 n e i is more important than K ee 
at all displayed densities for T = 10 s and 10 9 K. In these 
and we do not show K e to simplify the 
figure. However, n ee dominates in the pure ground-state 
matter for T = 10 7 K at p > 10 12 g cm -3 (in the inner 
crust of a neutron star). This dominance is fully pro- 
duced by the Landau damping of transverse plasmons in 
the ee-collisions. For p > 10 13 g cm -3 in this case we 
have K e « K ee . For an impure cold-catalyzed matter n e i 
remains important for all densities at T = 10 K, but 
K e e is comparable to K e i in the inner neutron star crust. 
Note that the values of K e i in Figure |3| take into account, 
in an approximate manner, the freezing of Umklapp pro- 
cesses in electron-phonon scattering at low temperatures 
(see, e.g., Ref. 6]). For T > 10 s K the freezing is unim- 
portant, but at T ~ 10 K it enhances K e j. A more rig- 
orous treatment of this freezing can partly remove this 
enhancement (A. I. Chugunov, private communication), 
which can somewhat decrease K e i and reduce the impor- 
tance Of Kee- 

Therefore, a correct treatment of ee collisions can con- 
siderably reduce the electron thermal conductivity in a 
cold neutron star crust, with T ~ 10 7 K. This can in- 
crease the time of heat diffusion from the inner crust to 
the surface |jg|. In particular, the effect can be impor- 
tant for the propagation of thermal waves produced by 
pulsar glitches and for the emergence of these waves on 
the pulsar surface. The emergence can be, in principle, 
observed and give valuable information on the nature of 
pulsar glitches [8|. 

While calculating n ee we have taken into account only 
the electron contribution into the polarization functions 
III and II t in Eq. © and neglected the ion contribution. 
As shown by Lampe |l3j | this is a good approximation 
at high temperatures T > Z 2 e 2 /(afce) at which the ions 



are weakly coupled and constitute a Boltzmann gas. The 
calculation of the ion polarization functions at lower T 
(in the regime of strong ion coupling, in an ion liquid 
or solid) is a complicated problem whose detailed solu- 
tion is still absent. We have also neglected the effects of 
strong magnetic fields which can greatly modify electron 
heat transport in strongly magnetized neutron star en- 
velopes. These effects are numerous (e.g., Yakovlev and 
Kaminker '26\) and are subdivided into classical (owing 
to rapid electron Larmor rotation) and quantum ones 
(owing to the Landau quantization of electron motion in 
a magnetic field). In principle, a proper treatment of ee 
collisions in a strong magnetic field can be performed in 
the same framework of the dynamical plasma screening as 
used above. However, the problem becomes much more 
complicated because the electron polarization tensor in a 
magnetic field is anisotropic (depends on the relative ori- 
entations of the wavevector q and the magnetic field) and 
cannot be generally decomposed into purely longitudinal 
and transverse parts |2l| . The effects of ion polarization 
and strong magnetic fields are outside the scope of the 
present paper. 



V. CONCLUSIONS 

We have reconsidered the electron thermal conductiv- 
ity K ee of degenerate electrons produced owing to the 
ee collisions taking into account the Landau damping 
due to the exchange of transverse plasmons (following 
the calculations of kinetic properties of quark plasma by 
Heiselberg and Pethick [17|L The Landau damping has 
been neglected in all previous calculations of K ee - We 
have analyzed the four regimes of the ee collisions in the 
degenerate electron gas (Section IIIII Table HJ and ob- 
tained analytic expressions for K ee which accurately ap- 
proximate the results of numerical calculations of K ee in 
wide ranges of the temperature and density of the matter. 
These results can be applied to study thermal conduction 
of degenerate electrons in metals, semiconductors, in de- 
generate cores of evolved stars and white dwarfs, and in 
envelopes of neutron stars. 

Our main conclusions are the following. 

(1) The Landau damping strongly modifies K ee in a 
relativistic degenerate electron gas, at densities p > 
10 6 g cm" 3 , but it is also quite noticeable at lower p 
(for instance, at p = 10 5 g cm -3 , see Figure^J. The Lan- 
dau damping increases the ee collision rate and decreases 
K ee , increasing the contribution of the ee collisions into 
the total electron thermal conductivity K e . 

(2) The most dramatic effect of the Landau damping 
on K ee takes place at p ^> 10 6 g cm~ 3 for temperatures 
T much below the electron plasma temperature T pe (the 
regime IV, Table [|J. In this case n ee shows a non- Fermi- 
liquid behavior; it becomes temperature independent and 
is described by the asymptotic expression (|2*U1) . 

(3) The conductivity K ee becomes comparable to 
the electron thermal conductivity K e ; provided by the 
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electron-ion collisions (and gives thus a noticeable contri- 
bution into the total conductivity K e ) in a warm plasma 
of light (low-Z) ions at temperatures T a few times lower 
than T pe (Figure |2J) ■ These conditions are typical for 
degenerate cores of white dwarfs and giant stars where 
thermonuclear burning of light elements can occur (par- 
ticularly, in the vicinity of the carbon ignition curve). 

(4) The conductivity re ee dominates over n e i and de- 
termines the total electron thermal conductivity K e of 
the dense pure cold-catalyzed matter at T ~ 10 7 K and 
P ~ 10 12 § cm -3 (in the inner crust of a cold neutron 
star, Figure 13). At these T and p the conductivity K ee 
is important even for impure cold-catalyzed matter. It 
can affect the propagation of thermal waves, excited in 
the inner neutron star crust during pulsar glitches, to the 
pulsar surface. 

The electron conductivity K e operates also in neutron 
star cores, at p > 1.5 x 10 14 g cm~ 3 . This conductiv- 



ity should also be reconsidered taking into account the 
Landau damping. Our present results cannot be directly 
applied to this case because in the cores the electrons col- 
lide efficiently at least with degenerate electrons, muons, 
and protons, and these collisions deserve a special study. 
We intend to analyze them in a subsequent publication. 
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